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Abstract 

The high-precision data from hadronic r decays allows one to extract information 
on QCD condensates. Using the finalized ALEPH data, we obtain a more rigorous 
determination of the dimension 6 and 8 condensates for the V — A correlator. In 
particular, we find that the recent data fix a certain linear combination of these 
QCD condensates to a precision at the level of 0(2) %. Our approach relies on 
more general assumptions than alternative approaches based on finite energy sum 
rules. 
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1 Introduction 



The physics of hadronic r decays has been an important testing ground for QCD since 
long. The precise data from the LEP collaborations have allowed us to obtain invaluable 
information on fundamental and effective parameters of the theory. With the comple- 
tion of the final analysis by the ALEPH collaboration [1] it is timely to exploit these 
data also for a more precise extraction of QCD condensate parameters. 

In a previous letter [2] we had used earlier ALEPH data [3] in a functional method 
to extract within rather general assumptions the QCD condensate of dimension D = 6 
related to the V — A current. Here we report the results of a reanalysis using the 
more precise final data from ALEPH. Moreover, we found that the higher quality of 
the experimental data allows us also to obtain information on the condensates of higher 
dimension. In particular, we find that the data imply a strong correlation between 
D = 6 and D = 8 V — A condensates. 

We extract the condensates from a comparison of the time-like experimental data 
with the asymptotic space-like results from theory. The assumptions of our approach 
are quite general. The essential property of the exact correlator in the space-like region 
on which our approach relies is that its fall-off with the 4-momentum squared s is 
determined by at most two lowest- dimension terms in the operator product expansion 
(OPE). When we aim at a determination of the D = 6 condensate, we assume a fall-off 
like 1/s 3 within an error band that scales like 1/s 4 . In the analysis where we want to 
determine the D = 6 and D = 8 condensates, we correspondingly assume a fall-off like 
a sum of 1/s 3 and 1/s 4 terms within an error band described by a 1/s 5 term. These 
assumptions are more general than those of other approaches since they do not refer to 
the positive axis in the s-plane. Furthermore our assumptions are quite independent of 
perturbation theory or indeed of QCD itself. Similarly to the Weinberg sum rules, they 
only depend on the pattern of chiral symmetry breaking. Of course, as QCD is well 
established, we like to discuss our results in the language of QCD and the OPE. For 
this reason we also include the 0(a s ) correction to the 1/s 3 term. If it was available, we 
could have included the perturbative correction to the 1/s 4 term as well. However, as 
we will see, the inclusion of higher-order terms does not qualitatively affect our results. 

One should, however, keep in mind that it is not possible (not even in principle) to 
reconstruct the correlator in the space-like region from error afflicted time-like data, as 
this constitutes an analytic continuation from a finite domain. One has to stabilize this 
"ill-posed problem" by suitable additional assumptions. For example, the application 
of finite energy sum rules (FESR) for an extraction of QCD condensates could be 
justified if the result of QCD and the operator product expansion in the space-like 
region is simply a series in powers of 1/s times condensates (vacuum expectation values 
of operators) and if there are no truly non-perturbative terms. In this case, each 
moment would pick out a single operator. However, this is unfortunately not the 
case since logarithms arising within perturbation theory do not fall in this class of 
functions. The higher and lower dimensional condensates contribute to a given moment 
starting with the inclusion of corrections of order 0(a 2 s ) [4]. It is also known that the 
perturbation series starts to diverge at some, not very high, order. As higher orders of 
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the relevant Wilson coefficients are not known, one can only hope that the extraction 
of the condensates is stable when correction terms are included. There is the additional 
problem of the contribution of the truly non-perturbative terms (non-OPE) to the 
integral on the circle in the complex plane, even if this uncertainty, which is expected 
to be most prominent near the physical cut, can be reduced by choosing suitable linear 
combinations of moments. 

We should expect that there is a price to be paid for the generality of assumptions 
of our approach. However, even if it turns out that the errors in the values of the 
extracted parameters are relatively large, we can hope that our results lend additional 
confidence to the numerical results obtained with the help of FESR. 

2 QCD condensates 

We consider the polarization operator of hadronic vector and axial-vector charged cur- 
rents, J ll = V fl = u^^d and J^ = A^ = uj^d, 

njL = i J dxe i(lx (TJ,(x)J u (Oy) (1) 

= (- M 2 + ?„<&,) n!rV) + ^<MV) . 

The conservation of the vector current implies Hy = 0. The connection to experimental 
observables is most easily expressed with the help of the spectral functions which are 
related to the absorptive part of the correlators. Keeping the normalization as defined 
in most of the experimental publications, the functions 

Vj{s) = 27rlmn{?(s), dj{s) = 2vrlmn^ } (s) (2) 

can be extracted from the decay spectrum of hadronic r-decays. We restrict the present 
study to the V — A component which is related to the branching ratios of r decays, 

B(t — > v T + hadrons, V — A) 
-tt-r.v-A = ^7 : ! , (3) 



through 



" l(s) -° l(s) = ^fcldTl 1 ^) [ 1 + 2 < ] • (4) 

ml di?A A s 



«o(«) = - — £ ^-r 1 T I (5) 

6\V ud \ 2 S EW B e ds\ ml) 

Here, B e = 0.17810 ± 0.00039 pQ is the branching fraction of the electron channel, 
is the CKM-matrix element, \V ud \ = 0.9746 ± 0.0006 [5], the r mass is m T = 1.777 
GeV and Sew = 1-0198 ±0.0006 accounts for electroweak radiative corrections [6]. The 
spin-0 axial vector contribution a Q (s) is dominated by the one-pion state, a Qi7T (s) = 
2n 2 f*5(s - ml), with the vr-decay constant f w = 0.1307 GeV [7]. 
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We use the final experimental data from the ALEPH collaboration [T] because they 
have the smallest experimental errorsS They are given by binned and normalized event 
numbers related to the differential distribution dR T y,Al 'ds. The quality of these data 
has increased a lot as compared to those of Ref. [3] which we used in our previous 
analysis |2]: first, the higher statistics has allowed the ALEPH collaboration to double 
the number of bins, and, secondly, the experimental errors decreased due to both higher 
statistics and a better understanding of systematic errors. 

The (V — A) correlator is special since it vanishes identically in the chiral limit 
(m q = 0) to all orders in QCD perturbation theory. Renormalon ambiguities are thus 
avoided. We will neglect in our analysis perturbative contributions proportional to the 
light quark masses. Non-perturbative terms can be calculated for large \s\ by making 
use of the operator product expansion of QCD 

n vSM = £ & i 1 + cS L0 ^ + ok 2 )) (6) 

D>6 V / \ / 

where O s [ ) ~ A are vacuum matrix elements of local operators of dimension D. For the 
correlator with spin + 1 and in the chiral limit, the sum begins with terms of dimension 
D = 6. Assuming vacuum dominance or the factorization approximation which holds, 
e.g., in the large- Nc limit, the matrix elements can be written as 



c 6 + In ( — 

— s 



8 V " A = -47ra s i(l-Na 2 )(qq)(qr P G^q); (7) 
O w w A = -lna s (qq) 2 (^m* + 32na a (G af) G'#) 



9 WJ/ V 9 

The complete expression for 0§~ A is known to involve two operators [!(]]. However, our 
analysis does not rely on the factorization hypothesis since we are going to determine 
Oq~ a and Og~ A , but not the condensates (qq) or (q^apG^q) separately. In the result 
for OYq~ A [H], ml is defined through the 5-dimensional quark-gluon mixed condensate. 
Starting from the second order, coefficients in perturbation theory depend on the reg- 
ularisation scheme implying that the values of the condensates are scheme-dependent 
quantities. The NLO corrections for 0^~ A were computed first in [12] and the coeffi- 
cient Cq was found equal to 247/12. This calculation was based on the BM definition of 
7 5 in dimensional regularisation. A different treatment of 7 5 as used in Ref. [X3J leads 
to c 6 = 89/12. 

The typical scales determining the condensates are around 300 MeV. For example, 
from the Gell-Mann-Oakes-Renner relation [13] one has (qq) ~ —(250 MeV) 3 , and 
from charmonium sum rules one expects {a s /n)(G 2 ) ~ (300 MeV) 4 [T3]. Taking the 



'Earlier measurements have been reported in [S], see also [5] and references given there. 
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results from FESR approaches as a guide, we expect also Og~ A and O^ ~ A to be of 
order 10 -3 GeV 8 and 10~ 3 GeV 10 , resp. (see e.g. [IB])- This is small enough so that 
the OPE makes sense. If O^ A would be much larger, radiative corrections to higher- 
dimension condensates would mix significantly with the lower-dimension ones through 
their imaginary parts. 

In the following we shall summarize the functional method^ introduced in Ref. [2] 
which allows, in principle, to extract the condensates from a comparison of the data 
with the asymptotic space-like QCD results under rather general assumptions. 

We consider a set of functions F(s) (to be identified with Ily^ (s)) expressed in 
terms of some squared energy variable s which are admissible as a representation of the 
true amplitude if 

i) F(s) is a real analytic function in the complex s-plane cut along the time-like 
interval Tr = [s ,oo). The value of the threshold sq depends on the specific 
physical application (s = (2m n ) 2 for U v , s = m% for IT4). 

ii) The asymptotic behavior of F(s) is restricted by fixing the number of subtractions 
in the dispersion relation between F(s) and its imaginary part along the cut 
f(s) = ImF(s + iO)\ s( zr R (for Ilyl^ (s) no subtractions are needed): 

F(s) = - ri^-dx + S., (8) 
■K J x- s s 

where the term with f% accounts for the contribution from the pion pole, (Eq. 
EJ), so that f(x) can be taken directly from the ALEPH data. 

We have two sources of information which will be used to determine F(s) and f(s). 
First, there are experimental data in a time-like interval T exp = [sq, s max ] with sq > 
for the imaginary part of the amplitude. Although these data are given on a sequence 
of adjacent bins, we describe them by a function / exp (s). We assert that f exp is a real, 
not necessarily continuous function. The experimental precision of the data is described 
by a covariance matrix V(s, s'). 

On the other hand, we have a theoretical model, in fact QCD. From perturbative 
QCD we can obtain a prediction for the amplitude in a space-like interval T L = [s 2 , si]. 
This model amplitude is a continuous function of real type, but does not necessarily 
conform to the analyticity property i). Since perturbative QCD is expected to be 
reliable for large energies, we expect that there is also useful information about the 
imaginary part of the amplitude provided that \s\ is large, i.e. we can also use QCD 
predictions for f(s) = lmF(s + «0)| se ( SmaxiOO ). To implement this idea, we rewrite the 
dispersion relation dSJ) in the following way: 

s 7T /, x — S 7T In x — s 

TThe functional method underlying our analysis has first been described in [T71 ITS] . 
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Then we can calculate the left-hand side of this equation in QCD, 



F QCD ( S ) := F QCD (s) - f l-- ^ ^Mrfx , (10) 

S TT L X — S 

u o max 

whereas the right-hand side can be determined from experimental data. Thus we shall 
test the hypothesis whether the left-hand side of (EJ) can be described by QCD while 
the right-hand side describes the experimental data. 

In order to do that, we need an a-priori estimate of the accuracy of the QCD 
predictions. This can be described by a continuous, strictly positive function ol(s) for 
s 6 Ti which should encode errors due to the truncation of the perturbative series and 
the operator product expansion. It is expected to decrease as \s\ — > oo and diverge for 
s — > 0. In Ref. [2] we aimed at a determination of the D = 6 condensate of IK/- a- 
There we were allowed to consider the contribution of the dimension D = 8 operator 
as an error, using o L (s) = \Oj~ A \ max /s 4: with |C^~ A | max in the order of 10~ 3 GeV 8 . If 
the perturbative part dominates, as is the case for the individual vector or axial vector 
correlators, the last known term of the perturbation series could be used as a sensible 
estimate of the error corridor, possibly combined with the first omitted term in the 
series over condensates. In the present work, we focus on obtaining information on the 
correlation of the condensates 0^~ A and 0^~ A ; therefore we will use in a similar way 
an estimate of |0}o~ A |max to define an error corridor which, in this case, scales with 
1/s 5 . 

The goal is to check whether there exists any function F(s) with the above analyt- 
icity properties, the true amplitude, which is in accord with both the data in r exp and 
the QCD model in I\. In order to quantify the agreement we will define functionals 
XlU] an d XrII] using an L 2 norm. For the time-like interval we simply compare the 
true amplitude f(s) with the data and use the covariance matrix of the experimental 
data as a weight function: 

X 2 R [f] = jp— | / dx / dx'V-\x,x') (f(x) - / cxp (x)) (/Or') - / cxp (x')) . (11) 

1 1 exp | J so J so 

Experimental data correspond to cross sections measured in bins of s, so that we can 
calculate this integral in terms of a sum over data points. The ALEPH data which we 
use are given for 140 equal-sized bins of width As = 0.025 GeV 2 between and 3.5 GeV 2 . 
x\ given in (jllj) is in fact the conventional definition of a \ 2 normalized to the number of 
degrees of freedom and has a probabilistic interpretation: for uncorrelated data obeying 
a Gaussian distribution we would expect to obtain \\ = 1. Since experimental data at 
different energies are correlated, we instead expect 

x-p = ^ E \lv(*i>*i)y(*i>*s)v~\*> s i) • ( 12 ) 

In order to define a measure for the agreement of the true function f(s) with theory, 
we use the left-hand side of flH]) which is well-defined and expected to be a reliable 
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prediction of QCD in the space-like interval for not too small \s\. This expression can 
be compared with the corresponding integral over the true function. Thus we define 

Xl[f] = pfrr / w L (x) (f qcd (x) - - f ^ ^-dx') 2 dx (13) 
Jr L V nJ S0 x'-x J 

where wl is the weight function for the space-like interval and identified with 1/<j\(s). 

In order to find a best estimate of the function f(s), we minimize xt subject to the 
condition xkif] — Xcxp- The solution can be obtained by solving a Fredholm equation of 
the second kind. This can be done numerically by expanding f(s) in terms of Legendre 
polynomials (see Ref. [2] for details). The algorithm to determine an acceptable value 
for the condensate is then the following: 

i) For the given value of xL P we determine the solution for f(s) iteratively and 
calculate the corresponding value of xlLf] as a function of 0^~ A . 

ii) We minimize this xtif] with respect to the values of 0^~ A . 

Hi) We determine the error on the condensates by solving x\ = Xl mm + ^X 2 where 
A% 2 is chosen in the conventional way to reflect 1 or 2a confidence regions for the 
values of C^ A . 

The size of a^ determines the minimal value of x\ according to Eq. (jT3"j) that can be 
reached by this algorithm. Obviously, widening the error corridor (increasing cr^) will 
lead to values for xi mm as small as desired. In such a case, the information obtained 
from the fit is not conclusive, since any model function f(s) can be made consistent with 
the data if one allows for a wide enough error corridor. On the other hand, narrowing 
the error corridor will increase xt mim signalling a bad fit, i.e. bad consistency of theory 
with data if the model function is not perfectly describing the data. However, a non- 
trivial result of our approach is the fact, that there exists a choice for the error corridor 
that leads to values ximm = 0(1). We shall assume that the underlying probability 
distribution is Gaussian, and choose xt mm an d accordingly <jl such that the fit result 
corresponds to a la CL. In practice this is done by adjusting \Oj~ A \ max (for the 1- 
parameter fit of 0^ A ) or |C^ A | max (for the 2-parameter fit of 0^ A and <D$~ A ). 



3 Numerical results and discussion 

We start with a discussion of results obtained by 1-parameter fits of the dimension 
D = 6 condensate. We find a consistent description of the data by QCD predictions 
including a non-zero value of Oq~ a as in our previous publication. The quality of the 
fit is good, i.e. a value of xt mm corresponding to a la CL can be fixed by choosing an 
error corridor described by the contribution from an 0^~ A condensate of the expected 
size of ~ 1.3- GeV 8 . A direct comparison of the data with the regularised function 
f(s) (see Fig. [T]) shows a nice agreement over the full range of s with the exception 
of the highest s-bins, where experimental errors are large. In comparison with [2], the 
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agreement between data and theory has improved. A number of consistency checks, as 
described in [2], has been performed and did not reveal any problem. 

The results of the 1-parameter fit can be summarized by quoting values for 0^~ A 
at LO and at NLO for the two available values of the NLO coefficients^: 

C%~ A = -5.9t\i x l(T 3 GeV 6 for £ 6 = 0, LO , (14) 
0^ A = — 4.9±i;f x l(T 3 GeV 6 for £ 6 = ^ , (15) 

Oj~ A = -3.6±i:l x 10- 3 GeV 6 for c 6 = ^ . (16) 

The NLO results are based on the 4-loop expression for a s with Ay^(Nf = 3) = 0.326 
GeV and the renormalization scale chosen equal to \s\. They are not sensitive to 
changing A^g within the present experimental error of ±0.030 GeV. Moreover, the fit 
results based on the two different values of the NLO coefficients c$ agree within errors 
and their difference with respect to the LO result is consistent with a shift calculated 
from the correction term choosing a typical value of 0(1.5) GeV for the renormalization 
scale in a s . This is to be expected since our method would work for any s-dependent 
ansatz for 0§~ A as well. 

The main result of our present analysis is obtained from a 2-parameter fit of the 
dimension D = 6 and D = 8 condensates. Here we do not include the NLO contribution 
to 0^~ A since the corresponding NLO coefficient for 0^~ A is not known. In this case we 
find agreement of theory and data at the la CL if we choose an error corridor described 
by the dimension D = 10 contribution with the value ~ 5.7 • 10~ 3 GeV 10 . The result 
presented in Fig. [2] shows a strong negative correlation of 0^~ A and 0^~ A . Both the 



*Note that the normalization differs by a factor of 2 with respect to [2] 
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Figure 2: 1- and 2a CL regions in the (O^ A , (9^ A )-plane from the 2-parameter fit 
(LO). The central values are marked by the dashed lines. 



central values as well as the errors from the 2-parameter fit are consistent with those 
from the LO 1-parameter analysis: The la range allowed for <D^~ A for fixed 0^ A = 
(which is the assumption underlying the 1-parameter fit) agrees with (1141) . However, 
leaving the value of 0^~ A unconstrained, as in the 2-parameter fit, one finds a much 
larger range for <D^~ A . The minimum value of x\ f° r the 2-parameter fit is located at 
the values 

0V-A = -6.8^| x 10' 3 GeV 6 , 

£>v- A = 2Atfl 7 x l(T 3 GeV 8 . (17) 

The errors on 0^~ A for fixed 0^~ A are small, but the allowed range for Oj~ A is not 
very restrictive (note the different scales for the two condensates in Fig. [2]). However, 
the strong correlation allows one to determine a linear combination of 0^~ A and Og~ A 
with a rather small error: 

Oj" A + 3.04 GeV 2 • 0^~ A = -18.307_°;| x lO^GeV 8 . (18) 

We also checked that we obtain consistent results when including the NLO correction 
to 0^~ A : The 1- and 2a contours are shifted, essentially without changing their form, 
to larger values of 0^~ A exactly as can be inferred from the 1-parameter fits: for cq = 
(89/12, 247/12) we find the minimum of xl at C^~ A x 10 3 /GeV 6 = -6.8 (-6.6, -5.8). 

There exist a number of previous extractions of QCD condensates in the literature, 
all based on a sum rule approach and at LO [H1H3I2Q1I2HI22I23I21I23I2B]- Th e 
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results for 0^~ A cover values between (-2.27±0.51) ■ 10 3 GeV 6 [22] and (-8.0±2.0) ■ 10 3 
GeV 3 [23] ■ These values correspond to a scale of about 200 MeV which is comparable 
to Aqcd- Although errors given by the authors are typically in the order of 25 %, their 
central values are only in rough agreement. The observed variations of these results 
represent the ambiguities inherent in the QCD sum rule approach. Our result nicely 
falls into the same range, also with an error estimate of the same size. 

For Oj" A , previous results range from (-10.8 ± 6.6) • 10~ 3 GeV 8 [2E] to (12.2 ± 
2.9) ■ 10" 3 GeV 8 [25]. A recent conservative estimate [21] is 0^~ A = (-I2± 7 n ) • 10~ 3 
GeV 8 . Again, our result agrees within the estimated precision. Errors for the D — 10 
condensate, O s [q A ) are typically larger and the spread of values found in the literature 
even larger: they range from (-18.2±5.9) • 10~ 3 GeV 10 [25] to (78±24)-10~ 3 GeV 10 [21]. 
These values are consistent with the value we had to choose for the error corridor to 
determine the la CL range in the 2-parameter analysis. Moreover, it is also interesting 
to note the agreement of the correlation between 0^~ A and O xr ~ A with corresponding 
results from [HI [25] . 

Obviously, our numerical results are, from a practical point of view, not superior to 
approaches based on finite energy sum rules. However, the fact that we find agreement 
within errors is not trivial. Since our approach is based on much more general assump- 
tions, the results obtained in this analysis give additional confidence to the numerical 
values obtained with the help of QCD sum rules. 
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